Charged Vortices in High Temperature Superconductors Probed by NMR 
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We report a first experimental evidence that a vortex in the high temperature superconduc- 
tors (HTSC) traps a finite electric charge from the high resolution measurements of the nuclear 
quadrupole frequencies. In slightly overdoped YBa2Cu307 the vortex is negatively charged by trap- 
ping electrons, while in underdoped YBa2Cu408 it is positively charged by expelling electrons. 
The sign of the trapped charge is opposite to the sign predicted by the conventional BCS theory. 
Moreover, in both materials, the deviation of the magnitude of the charge from the theory is also 
significant. These unexpected features can be attributed to the novel electronic structure of the 
vortex in HTSC. 



I. INTRODUCTION 

One of the most important physical properties of the 
vortex created in the type-II superconductors is that a 
vortex line can support a magnetic flux with a flux quan- 
tum $ = hc/2e(= 2.07 x 10~ 7 Oe-cm 2 )El. This fact has 
been confirmed experimentally more than 4 decades ago. 
On the other hand, it is only very recent that another 
prominent feature, namely the possibility that a vortex 
of the superconductor can accumulate |-a.-|finite electric 
charge as well, has come to be realizedaEl. The vortex 
charge appears as a result of the chemical potential dif- 
ference between the vortex core and the region away from 
the vortex core. It should be emphasized that the sign 
and magnitude of this charge is closely related to the mi- 
croscopic electronic structure of the vortex, which in turn 
reflects all the fundamental natures of the superfluid elec- 
trons and the low energy excitation out of the condensate. 
Moreover, it has also been pointed out recently that the 
vortex charge, strongly affects the dynamical properties 
of the vortextThj. For instance, the origin of the vortex 
Hall anomaly has been attributed to the vortex charge. 
Though the clarification of the issue of the vortex charge 
serves as an important test of the predictions for the vor- 
tex electronic structure and the dynamics, it has never 
been examined experimentally so far. This is mainly be- 
cause in conventional superconductors the magnitude of 
the accumulated charge within the core is very small and 
is extremely difficult to observe. 

In this paper we report a first straightforward attempt 
to identify the vortex charge in the high temperature su- 
perconductors (HTSC) by high resolution measurements 
of the nuclear quadrupole frequency vq, which is very 
sensitive to the local charge density. We show that a 
vortex in HTSC indeed traps a small but finite elec- 
tronic charge as well. In slightly overdoped YBa2Cu307 
the vortex is negatively charged, while in underdoped 
YBa2Cu40g it is positively charged. The sign of the 
trapped charge is opposite to the sign predicted by the 
conventional BCS theory. Moreover, in both materials, 



the accumulated charge is much larger than expected in 
the ordinary superconductors. We discuss several possi- 
ble origins for these discrepancies. 

II. VORTEX CHARGE 

We first briefly introduce the vortex charge in type-II 
superconductors. In the core of the conventional super- 
conductors, the distance between each discrete energy 
levels of the quasiparticles in the Andreev bound states 
are merely of the order of a few mK. Thus it is sufficient to 
view the energy levels as forming a continuous spectrum, 
just like in a normal metallic statetlO. Generally the 
chemical potential /i in the superconducting state differs 
from that in the normal state, if an electron-hole asym- 
metry is present. Assuming therefore that the vortex core 
is a normal metallic region surrounded by the supercon- 
ducting materials, this difference in \x is expected to arise 
and should lead to the redistribution of the electrons. In 
order to maintain the same electrochemical potential on 
both side, the charge transfer occurs between the core 
and the outside. 

In the framework of the BCS theory taking into ac- 
count the metallic screening effect, the charge accumu- 
lated within the vortex core per layer normal to the 
magnetic field is given as, 

where Xtf is the Thomas- Fermi screening length, s is the 
interlayer distance, fixs the chemical potential and e(> 0) 
is the electron chargeH. The sign of the core charge is de- 
termined by dlnT c /c?ln/x, which represents the electron- 
hole asymmetry. Outside the core, the charges with op- 
posite sign screen the core charge, similar to a charged 
particle in a metal. Far outside the core, this screening 
charges decay gradually with a power law dependence as 
r- 4 (seeFig.l (a))il. In strong fields (H cl < H <C H c2 ), 
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FIG. 1. (a):Schematic figure of the charge distribution 
around a single vortex core when the electrons are trapped 
within the core (negatively charged core). A(r) is the super- 
conducting energy gap. p(r) is the charge density. The charge 
accumulated inside the core are screened by the charges with 
opposite sign. p(r) decays gradually as r -4 well outside the 
core region. Q(r) is the total charge within the distance r. 
Q(r) goes to zero as r — » oo due to the requirement of the 
overall charge neutrality, (b): The charge density modula- 
tion in the strong magnetic field (H c i <C H <C H C 2) where 
each vortices overlaps. 

each vortex overlaps with its neighborhood. Then the 
charge density outside the core is nearly constant, and 
the periodic modulation of the charge density appears 
for the periodic vortex lattice (Fig. 1(b)). In ordinary 
superconductors, |Qf | is estimated to be ~ 10~ 5 — 10 _6 e, 
using k F - lA" 1 , X TF ~ hp 1 ~ lA, £ ~ 100A and 
\dlnT c /dliifi\ w \n(hoj£)/kBT c ) ~ 1 — 10, where ujd is 
the Debye frequency. Thus \Q$\ is negligibly small and 
is very difficult to observe. 

However, the situation in the case of HTSC seems 
to be promising because £ is extremely short compared 
to that of the conventional superconductors. More- 
over the strong electron correlation effects and d-wave 
pairing symmetry of HTSC are expected to chaiLgfi the 
electronic structure of the vortex dramaticallyLjllZl. In 
fact, recent STM measurements revealed that the vor- 
tex of HTSC is yery-different from those of conventional 
superconductorsE3~E2l. These unusual features of HTSC 
are expected to enhance the charging. We will discuss 
these issue later. 

The vortex charge also plays an important role for the 
vortex dynamics. When a vortex moves in the super- 
fluid electrons, the core plays a key role in dissipation 
processes!!!! One of the most striking phenomena is the 
vortex Hall anomaly, namely sign reversal of the flux flow 
Hall effect below T c , which is observed in most HTSCtil. 
This Hall sign reversal indicates that the vortices move 
upstream against the superfluid flow. Such an unusual 
motion has never been observed in any other fluid includ- 
ing superfluid helium and cannot be explained in the 
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FIG. 2. Doping dependences Q^-iVq on the Cu(2) site, for 
YBa 2 Cu 3 7 _5 Jma^Sr.CuO^ and HaBa 2 Cu0 4+a .Ea v Q 
is proportional to the hole number. 

framework of the classical hydrodynamic theoryS. Re- 
cently, this phenomena has been discussed in terms of the 
vortex charge which produces an additional force acting 
on the vorticesoia. 

To elucidate the vortex charge, a direct observation of 
the local carrier density is strongly required. It has been 
pointed out by many authors that-^a-in HTSC is very 
sensitive to the local hole densitynJO. In a solid, an 
electron distribution with spherical asymmetry such as 
unclosed 3d shell and noncubic surrounding ions induce 
a local electric field gradient (EFG) in the vicinity of the 
nuclei. This EFG lifts the degeneracy of the nuclear spin 
levels, interacting with the nuclear quadrupole moment, 
Qn- The relevant information is obtained from the nu- 
clear quadrupole resonance (NQR) in zero magnetic field 
and the nuclear magnetic resonance (NMR) in a finite 
magnetic field. For the 63 Cu nuclei with spin 1=3/2, the 
NQR resonance frequency Vq® R is expressed as 
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where eq zz is the largest principle (z— axis) component 
of EFG at the nuclear site, Q N (=-0.211 barn for 63 Cu) 
is the quadrupole moment of copper nuclei. Ell The asym- 
metry parameter rj of the EFG defined as 77 = \(q xx — 
q yy )/q Z z\ is close to zero for the Cu site in the two di- 
mensional Cu02 planes (Cu(2) site). 

In a strong magnetic field when the Zeeman energy is 
much larger than the quadrupole energy, each Zeeman 
level is shifted by the quadrupole interactions and thus, 
two satellite peaks (±3/2 <-> ±1/2) appear on both sides 
of the central (±1/2 <-> Tl/2) resonance peak. The fre- 
quency difference between the upper and the lower satel- 
lites exactly coincides with 2vq, when a magnetic field 
is applied parallel to the largest principle axis of EFG, 
namely, H \\ c-axis in the present experimental 
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FIG. 3. The NQR (a), the lower (b) and the upper (c) 
NMR satellite spectra at 9.4T for 63 Cu(2) site of YBa 2 Cu 3 7 
for various temperatures. (220K, 200K, 180K and 160k to 10K 
by a 10K step.) 

condition described below. Even for the second order or 
more higher order quadrupole effect in the presence of 
the asymmetry r\ of EFG, the each satellite line shifted 
as much and the corrections for 2vq are vanishing in the 
case of H || c-axis.cJ. 

Generally the EFG originates from two different 
sources, namely from the on-site distributions q on - S ite 
of the electrons and from the surrounding ions qi on , 
q = Qon-site + qion - Recent analysis of q on the Cu(2) 
site suggests that (j „_ S ij_-.is mainly composed of the 
Cu 4p and 3d shell termsE3. In HTSC the holes in the 
Cu 3d x 2_ y 2 orbital play an important role for the on- 
set of superconductivity. Figure 2 shows the doping de- 
pendence of ^o-pf the Cu(2) site for._yBa2Cu307_,5Ea, 
La 2 _ K Sr :c Cu04E3, and HaBa 2 Cu04 + 5.E£l In all materials, 
vq increases linearly with the number of holes in the 
planes and can be written as, 

v Q = An hoie + C, (3) 

where rihoie is the number ©Wholes per Cu(2) atom, and 
A and C are constantsE-aES. Although C is strongly 
material dependent, reflecting the difference in 2/j on , A 
w 20 — 30 MHz per hole for the Cu(2) atom is essentially 
material independent. Thus, the precise measurement 
of vq makes possible an accurate determination of the 
change of the local hole number at the Cu site. 

III. EXPERIMENTAL 

The principle of our experiment is the following. In 
the measurement, only the resonance of the 63 Cu(2) nu- 
clei outside the vortex core is detected. This is because 
the applied field is much less than H C 2 and hence the 
core region occupies a smaller area in the sample. If the 
vortex core traps (expels) a finite amount of electrons, 
the electron 
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FIG. 4. The NQR in zero field (a) and the lower (b) and 
the upper (c) NMR-satellite spectra at 9.4T for 63 Cu(2) site 
of YBa 2 Cu 4 8 for various temperatures. (300K, 220K, 200K, 
180K and 160k to 10K by a 10K step.) 

density outside the core should decrease (increase) from 
that in zero field where the electron distribution is uni- 
form, as shown in Fig. 1(b). We are able to detect the 
change of carrier density through the change of the value 

of VQ. 

In the present measurements, we used slightly over- 
doped YBa2Cii307 and underdoped YBa 2 Cu40s in 
which the NQR and NMR spectra are very sharp com- 
pared to those of other HTSC. The NMR spectra are 
obtained for fine powdered sample (the grains are less 
than 33 /im) with uni-axial alignment. The each grain 
aligns to an easy axis (the c-axis) in a high external field 
at room temperature, and then fixed with epoxy (Stycast 
1266). 

The Cu-NQR and NMR spectra were obtained by a 
conventional pulse spectrometer. The NMR experiments 
were performed in the field cooling condition under a 
constant field of 9.4T by using a highly-homogeneous su- 
perconducting magnet which was stabilized to less than 
lppm during the experiment. The present measurement 
in the vortex state is made in the so-called Bragg-glass 
phase in which the quasi-long-range order of the vortex 
lattice is preserved. We see slightly-asymmetric shape 
of the NMR-satellite spectra due to the somehow miss- 
alignment of some grains. The line broadenings of the 
NQR spectra in zero external field are not observed even 
below T c in both YBa 2 Cu 3 7 and YBa 2 Cu 4 8 . The 
width of the NMR spectra becomes broader with decreas- 
ing T below T c partly due to the magnetic field inhomo- 
geneity caused by the introduction of the vortices and 
partly due to the spatial distribution of carrier density, 
which we will discuss later. 
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FIG. 5. The temperature dependence of vq obtained from FIG. 6. The temperature dependence of vq obtained from 
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IV. RESULTS AND DISCUSSION 

Figure 3 shows the NQR and lower and upper 
satellite spectra of 63 Cu(2) NMR for optimum doped 
YBa 2 Cu 3 7 . Figure 4 shows also the NQR and NMR 
spectra of 63 Cu(2) site of YBa 2 Cu40s. All data were 
taken by utilizing a phase coherent spin echo spectrome- 
ter. The spectra are obtained with a superposed method 
of the Fourier Transform spectra of the spin echo mea- 
sured at a certain frequency interval. The temperature 
dependences of vq obtained fron NMR and NQR are 
shown for YBa2Cii307 and YBa2Cu40s in Figs. 5 and 6, 
respectively. The temperature dependences of Vq® r for 
both YBa2Cu3C>7 and YBa2Cu4Q« are quite similar to 
those of fq reported previously^. It should be noted 
that the procedure for obtaining vq from the frequency 
difference between the upper and the lower satellites is es- 
sentially free from the influence of the change of magnetic 
shift (or Knight shift). Moreover, the magnetic effect of 
the asymmetric broadening (so-called Redfield pattern) 
due to the vortex lattice in the superconducting state is 
exactly cancelled out in the process determining vq from 
NMR. Thus, we obtained the vq values simply from the 
difference of the peak frequencies of the two satellite lines. 

We plot in Fig. 7 the difference between vq in zero 
field and in the vortex state, Avq = vq(Q) — vq(H), 
for YBa2Cu3C>7 and YBa2Cu40s- The vq in zero field 
for NQR is o btained after the correction by the factor 
of \/l + 1] 2 /3 in Eq.(l), although this factor is at most 
0.03% of vq for r\ «0.04 of the present materials. In both 
materials Avq is essentially zero above T c , indicating no 
modulation of the carrier density. Meanwhile a 



nonvanishing Avq is clearly observed below T c in both 
materials. While Avq ~ —25kHz is of negative sign 
and in YBa2Cu3C>7, Avq ~ 50kHz is of positive sign 
in YBa 2 Cu 4 8 at T=0. 

We discuss here several possible origins for the nonzero 
Avq. We first point out that the magnetostriction can- 
not be the origin of the nonzero Avq. In fact, both 
the magnetostriction and ultrasonic absorption measure- 
ments showed that the local lattice distortion A£ caused 
by the magnetostriction is negligibly small under the field 
cooling condition; A£/£ < 10 -8 below 10 T where £ is 
the lattice constant. We next remark on an in-plane 
charge modulation caused by a charge stripe formationE2l 
or th&daarge density wave (CDW) transition of the chain 
site.EJt3 Although the static charge ordering associated 
with a stripe formation has been discussed in some of 
the high-T c cuprates, such orderings have never been re- 
ported in YBa2Cu 3 07 nor in YBa2Cu40s. Moreover, it 
has been reported that in La:214 compounds the static 
charger-stripe order gives rise to the wipe-out of NQR 
signalsES, implying that the static stripe order, if present, 
induces a large effects on quadrupole interactions. How- 
ever, we have not observed the broadening of the NQR 
spectra (in zero field) and nor any sign of the wipe-out of 
NQR and NMR signals in the present measurements for 
YBa 2 Cu 3 07 or YBa2Cu40s- Finally we mention a recent 
report of the broadening of the NQR line widths in the 
plane at low temperatures in YBa2Cu3 07_,5, which was 
explained in terms of the CDW formation in the chain 
site.E3 In the present study, in which we obtained sh arpe r 
NQR line widths compared to the previous reportsEJE3, 
such line-broadenings due to the electric quadrupole in- 
teraction is not observed. Thus there is no evidence of 
the charge modulation due to CDW or stripe formation 
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FIG. 7. T-dependence of Avq = u Q (0) - u Q (H = 9.4 T) 
for YBa2Cu307 and YBa2Cu40g. In both materials nonzero 
Avq is clearly observed below T c , showing that the electron 
density outside the core differs from that in zero field. 

in zero field in our crystals. It is also quite unlikely that 
the magnetic field induces the CDW or stripe formations. 
Having ruled out these various possibilities, we conclude 
that the nonzero values for Avq naturally lead to the 
fact that the electron density outside the vortex core is 
different from that in zero field. 

We now discuss the issue of the sign and magni- 
tude of the accumulated charge. The negative Avq in 
YBa2Cu307 indicates the increment of the hole density 
outside the core. This excess density of holes is nothing 
but the holes expelled from the core. Therefore the accu- 
mulated charge in the core of YBa2Cu30y is negative. By 
the same reasoning, the positive Avq in YBa2Cu40s in- 
dicates a positive accumulated charge. Meanwhile, since 
the chemical potential decreases monotonically with dop- 
ing holes, Eq.(l), sgnQj =sgn(dlnT c /dln/i), predicts 
that > in the underdoped regime while < 
in the overdoped regime. This is strikingly in contrast 
to the sign determined by the present experiment. The 
deviation of the magnitude of the charge from theory 
is also noteworthy. The magnitude of charges per pan- 
cake vortex, which are roughly estimated using ss 
Av Q H c2 /AH assuming H c2 ~ 200 T are Q c - -0.005e 
to -0.02e for YBa 2 Cu 3 7 and Q 6 ~0.01e to 0.05e for 
YBa2Cu4C>8. However, according to Eq.(l), is esti- 
mated to be ~ 10~ 4 — 10 _5 e where we assumed £ ~ 30A. 
Therefore \Qt\ determined by the present experiments 
are still one or two order of magnitude larger than ex- 



pected by Eq.(l). Thus the BCS theory not only predicts 
the wrong sign of the charge but also underestimates \Q$\ 
seriously. 

There are several intriguing possible origins for these 
discrepancies. For example, because of the extremely 
short £, the vortex core may be in the quantum limit 
kpt; ~1, where kp is the Fermi wave number. In this 
limit, the description of the quasiparticles in terms of 
semiclassical wavepackets bxeaks down in contrast to con- 
ventional superconductorsEl Furthermore, as suggested 
by recent theories of the vortex core based on e.g. the 
t-J or SO (5) models, the antiferromagnetic (AF) state 
may be energetically prtjfatable to the metallic state in 
the vortex core of HTSOi^liaO. If this is indeed so, the 
AF correlation is expected to enhance the charging effect 
because it causes a large shift of /j, by changing the den- 
sity of states of the electrons inside the core dramatically. 
We note here that the present results-exclude the possi- 
bility of the SO (5) insulating AF corerJO in which holes 
should be expelled from the core and the accumulated 
charges are always negative; the present result yields the 
opposite sign for YBa 2 Cu40s in the underdoped regime 
where the AF correlation is important. Therefore a de- 
tailed microscopic calculation is needed to evaluate the 
accumulated charge quantitatively including the sign. 



V. CONCLUSION 

From the precise 63 Cu-NMR and NQR measurements, 
we have shown that a vortex in type-II superconductors 
can trap a finite electric charge as well as magnetic flux. 
In the slightly overdoped YBa2Cu307 the vortex is neg- 
atively charged, while in underdoped YBa2Cu40s it is 
positively charged. In both high T c materials, the accu- 
mulated charges are much larger than expected in the 
ordinary superconductors. The sign and value of the 
charges indicate the novel electronic structure of the vor- 
tex in HTSC. 
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